Numerical study of induced vortex tunneling 
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Tunneling of vortex-antivortex pairs across a superconducting film can be controlled via inductive 
coupling of the film to an external circuit. We study this process numerically in a toroidal film 
(periòdic boundary conditions in both directions) by using the dual description of vortices, in which 
they are represented by a fundamental quantum field. We compare the results to those obtained in 
the instanton approach. 

PACS numbers: 85.25.Hv, 03.67.Lx 



m 
o 
o 

(N 

» 

o 
O 

m 
> 
o 
l> 
o 



(N ■ 

o : 
^ ■ 

^— > ■ 
c ■ 

s : 
cr 



X 



I. INTRODUCTION 

Persistent-current superconducting devices, in which 
the basis states are characterized by different vàlues of 
the enclosed flux, are interesting physical systems in their 
own right and are also promising candidates for applica- 
tions to quantum memory. The possibility to form quan- 
tum superpositions of macroscopic flux states has been 
demonstrated in experiments with SQUIDs [1,2]. 

Once reliable storage of quantum superpositions is 
achieved, it becomes necessary to consider possible mech- 
anisms for reading and writing operations and for assem- 
bling several such individual units (qubits) into a scalable 
quantum computer. In SQUIDs, various proposals have 
exploited the existence of a potential barrier between two 
stable basis states and have involved manipulation of the 
barrier itself, use of tunneling effects, or a controlled ex- 
citation over the barrier [3,4]. 

In this paper, we examine a model which describes tun- 
neling of vortices (short Abrikosov flux lines) in a ring 
of thin superconducting film. We consider a scheme in 
which a pulse of supercurrent suppresses superconduc- 
tivity, thus lowering the potential barrier and inducing 
tunneling, see Fig. 1. Because this process changes the 
flux enclosed by the ring, it can be used to form arbitrary 
superpositions of the basis states. We will see that in a 
suitable geometry, it is possible also to independently 
control the energy bias between the basis states, simi- 
larly to how it is done in SQUIDs. On the other hand, 
this device does not contain any Josephson junctions, 
thereby avoiding dissipation due to various fabrication 
issues, such as defects in the insulating barrier. 

Motivated by these observations, we consider a simpli- 
fied model, convenient for numerical work, in which the 
film has periòdic boundary conditions in both directions, 
forming a torus. The suppression of superconductivity 
and the biasing flux are represented in this model by two 
independent parameters: the vortex pair-production fre- 
queney M(t), and the driving force F(i). 

For applications to quantum computing, of main in- 
terest is the adiabatic regime, when there is very little 
residual excitation left in the system after the pulse (in 
other words, no real, as opposed to virtual, vortex pairs 




are produced). If this condition were not satisfied, the 
remaining vortices would be easily "detected" by the en- 
vironment (e.g. by electrons at the vortex cores), and 
that would result in ràpid decoherence. Thus, we envi- 
sion a situation when a virtual vortex and an antivortex 
are created, say, on the inside of the torus, transported 
along the opposite semicircles to the outside, and anni- 
hilated there, almost without a trace. 



L. 



L y 

FIG. 1. A schematic of vortex tunneling induced by a pulse 
of supercurrent I(t). The circles denote a vortex and an an- 
tivortex. 

In thin films, the elàstic mean-free time of electrons is 
short [5]; for a film of thickness d we use r e i = 2d/vp. 
This results in a strong suppression of the Magnus force 
on the vortex and a relatively small friction [6-9]. In 
addition, and in contrast to motion of real vortices, the 
density of normal electrons at the vortex core during tun- 
neling is a variational parameter, which adjusts itself to 
maximize the tunneling ratc. This lcads to further can- 
cellation of the Magnus force and a reduction in the in- 
ertial mass of the vortex. 

There are two theoretical approaches to induced vor- 
tex transport. One is based on instantons, solutions to 
the Euclidean equations of motion. Using the expres- 
sions [6,8] for the friction caused by the core fermions, 
we find that it results in the following contribution to 
the Euclidean action: 
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where n e = k F /3ir 2 is the electron density, L x is the 
width of the film, and TlüJq ~ A 2 /ep. Using kp = 1 À -1 , 
2lüq = 10~ s kpVF, L x = 1/im, and d = 4 nm, we obtain 
Sf = 170. By itself, this is not a small number, but the 
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crucial point is that Sf depends quadratically on the gap 
A. So, when we suppress A by a pulse of current, we 
also reduce Sf. In fact, the dependence of Sf on A is 
precisely the same as that of the vortex pair-production 
frcquency M. So, in what follows we simply includc the 
cffcct of the friction in our dcfinition of M. 

The driving force F, due to the energy bias between 
the basis states, can be viewed as a Lorentz force caused 
by an effective elèctric field, E', in the y direction. Thcn, 
the average vortex current due to tunncling is obtained 
in the instanton approach as [10] 
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E = (d/4e)E', and e the magnitude of the electron 
charge. In (2), the first term is due to instantons, and 
the second to anti- instantons; S'ò = ML x /c\ is the real 
part of the instanton action, and c\ is the limiting speed 
of vortex motion. 

The second method is entirely real-time. Vortices are 
described by a fundamental quantum field, and a nonzcro 
average vortex current comes out as a result of the dis- 
creteness of field modes. The periodicity in E has been 
confirmcd analytically in this approach [10], provided the 
vortex-antivortex potential can be replaced by its aver- 
age and included as an additional contribution to M. The 
second principal effect seen in (2) — the exp(— ML x /c\) 
dependence — has been confirmed only for the case of 
small E, E <C 2n/L x . It is of interest to develop this 
approach further, so that it can be applied also to cases 
with large E and a non-trivial potential. The present pa- 
per addresses, via numerical integrations, the first part 
of this program. 

In practice, it may be easier to fabricate a thin strip 
than a thin cylindcr. The field-theoretical method will 
apply to that case as well, provided one can establish the 
boundary conditions for the vortex field at the edges of 
the strip. 

The paper is organized as follows. In Scct. II we dis- 
cuss how one can independently control the vortex mass 
and the driving force in a thin superconducting ring. In 
Sect. III, we discuss the real-time description of vortex 
tunncling. Numerical results are presented in Sect. IV, 
and a summary in the concluding Sect. V. 



II. CONTROL OF CURRENT AND FLUX IN 
THIN SUPERCONDUCTORS 

Consider a uniform superconducting ring inductively 
coupled to an external circuit. Supposc the order param- 
eter winds n times around the ring. Then, the London 
current can be expressed through the flux $ enclosed by 
the ring as 



/ = -c($-n* )/* 



(3) 



Here $o is the flux quantum, and í is the "kinetic" 
ductance: 



mc L v 



tó 1 



(4) 



n s is the density of superconducting electrons, S is the 
cross-sectional area, and L y is the length of the ring. 

We have introduced the London penetration depth A, 
and because superconductivity can be deliberately sup- 
pressed by some means, A is in general different from the 
unperturbed penetration depth A. Recali also that in a 
thin film A determines the strength of the London cur- 
rent, but not the screening length of the magnètic field 
[11]. 

The flux $ in (3) is the total flux, which consists of the 
flux created by the external circuit and that created by 
the current / itself: 



$ = $c X t + -Col 
c 



(5) 



where C is the ordinary inductanec of the ring. Using 
this together with eq. (3), we can express the supercur- 
rent through $ e xt : 



n$ 



(6) 



Even though i and Cq enter eq. (6) symmctrically, there 
is an essential difference between them, since C$ depends 
only on the geometry and in this sense is a constant, 
while í depends on n s . 

So, there are two distinct limits of eq. (6). If i <C 
Co, i.e., the cross-sectional area S is large enough, then 
according to (6) <I> 0X t controls the current. On the other 
hand, if l » Cq, $ ox t determines only the product £1, 
i.e., the ratio I/n s . 

In the second, thin-ring, regime, by a slow variation 
of I/n s we can smoothly change the order parameter tp 
from a large initial value, for which vortex tunneling will 
be strongly suppressed, to some much smaller vàlues that 
allow tunneling, and then back to the initial state. For 
small, slowly changing ip, this can be seen directly from 
the Ginzburg-Landau (GL) equation: 



-H -i + r>l 2 = o 

en s v cr I \a\ 



(7) 



where j is the current density, v cr = h/2m£, is a critical 
velocity (£ is the coherence length), and a, b are GL pa- 
rameters. If quantum coherence can be preserved during 
this process, such a device would be reliable quantum 
memory. 

Now, £ ~ 2Ly\n.{L y /L x ), so the crossover between 
the thick and thin-ring cases oceurs at 



S cr ~27rA 2 In- 1 (L y /L x ) 



(8) 



Assuming that the logarithm is of order unity, and using 
an unperturbed value A = A = 150 nm, we obtain 



2 



1.4 x 10 5 nm 2 . 



(9) 



When n s (which is proportional to |V"| 2 ) is suppresscd to 
allow tunneling, A and, consequently, S CT grow, so if the 
condition S < S CT held in the initial state, it would hold 
even better during tunneling. 

According to this estimate, if the ring is made from a 
thin film, it does not have to be particularly narrow to 
achieve the thin-ring condition S < S C r- For examplc, for 
thickncss d = 4 nm, the estimate (9) allows for widths as 
large as 10 fim. 

As a specific example of how n s can be suppressed 
by a pulse of current, consider the double-arm geome- 
try shown in Fig. 1. It is convenient to imagine that 
the wire carrying the constant current / is closed at a 
large distance, so that the device can be viewed as a su- 
perposition of two closed circuits, with currents I\ and 
I = li + I<i. In addition to kinetic inductances £\, £2, 
the circuits have self-inductances £n, £22 and a mutual 
inductance —£12. We can also definc C\ = £n — £12 
and C2 = £12- Consider regime when the inductance of 
arm 1 is mostly kinetic, l\ ^> £1, while that of arm 2 
is mostly ordinary, l 2 <C £2- Then, the currents in the 
arms are 



h = 7 — [C2I - c($ cx t - n$ )] , 

■í-tot 

h = 7^-[^i/ + c($ext-n$ )] , 

i-tot 



(10) 

(11) 



where £ to t = í\ + C 2 - Suppose further that l\ 3> £2. I n 
this case, we see from (10) that the current I controls the 
parameter I\£\, which according to (7) determines how 
close the first arm is to criticality. 
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FIG. 2. A double-arm device, in which suppression of su- 
perconductivity in the weaker (first) arm is controlled by an 
external current /, while a biasing flux controls the energy 
diflerence between two flux states. 



The inductive energy of the double-arm device is 
1 



£ = 



2£ t 



-($ oxt -n$ r + e(n 



(12) 



where e(I 2 ) is independent of n and the external flux. If 
the device is biased by half a flux quantum, <í> ex t = $0/2, 
the energy (12) has two degenerate minima at n = 0, 1. 
If <& cxt deviates from half-quantum by a small amount 



A$ ext , the energy difference between the two minima is 
A£ = $ A"I> ox t/£tot· This results in an additional force 



F = A£/L x 



(13) 



acting on a vortex as it tunncls across arm 1. Notc that 
unlike the case of a single current-biased superconducting 
wire [12], this force is not rclatcd to the total current / 
but is an indcpendcntly controlled parameter. The only 
restriction is that |A£| should not exceed the energy 2M 
of vortex pair production, so that no real vortices are 
able to nucleate. 



III. REAL-TIME DESCRIPTION OF VORTEX 
TUNNELING 

Motivated by the arguments of the preceding sections 
we consider a model of vortex tunneling, in which the 
main force acting on the vortex is the driving force 
(13). The requisite suppression of the order parameter is 
achieved by some independent means, such as a pulse of 
external current in the double-arm device. 

In the real-time approach, vortices are described by a 
quantum field \ [13], which in our case obeys the equa- 
tion of motion 



d t 



-U(x,t) 



X 



(14) 



-c\dlx + M\t) X = Q. (15) 



Here d is the thickness of the film, e is the magnitude of 
the electron charge, and E(t) is a time dependent elèctric 
ficld produced by the changes in / and $ oxt . The speed 
ci is the limiting speed of vortex motion: in the second- 
quantized description (15) it plays the role analogous to 
the speed of light in relativistic quantum theory. 
The driving force is represented by the potential 



U(x,t) = -F(t)a 



(16) 



Such an explicit x-dependence in the equation is incon- 
venient for numerical work. However, it is possible to 
make a transformation of the field x, similar to a gauge 
transformation, so that the force disappears from the first 
term in (15) and appears instead as an addition to the 
elèctric field E: 



X 

u 



X exp[- J U(x,t')dt'} , 



4e /"* 

E -> E' = E - — / F(t')dt' . 



(17) 
(18) 
(19) 



In what follows, we will use the same notation \ f° r th° 
transformed field as we used for the original one. 

In general, the transformations (17)— (19) lead to one 
spurious effect. Imagine that F(t) starts from zero at 



3 



t = — oo, goes through nonzero vàlues near t = and thcn 
back to zero at t — oo. Then, according to eq. (15), U has 
no effect at t — > oo, while according to (19) the correction 
to E is still nonzero (and proportional to the integral of 
F). Fortunately, owing to the periòdic dependence of 
vortex transport on E, cf. eq. (2), this correction is 
immaterial provided F satisfies a quantization condition: 



J — i 



F(t)dt = 



2-Khrí 



(20) 



where n' is an integer. Only in this case the problem 
obtained by the transformations (17)— (19) is equivalent 
to the original problem (15). 

The transformed equation has no explícit dependence 
on x and can be solved by the mode expansion 



X (x,í) = VhY, [«k/k(í) +/?l k /k*(í) 



ikx 



(21) 



where k = (k x ,k y ), a and j3 are annihilation opcrators 
for vortices and antivortices, respectively, and /k(í) are 
the mode functions that take into account the time de- 
pendence of E'(t) and of the vortex "mass" M(t). We 
have assumed that the vortex ficld has periòdic boundary 
conditions in both dircctions. 

Substituting the expansion (21) into the field equation, 
we obtain the equation for the mode functions: 



/ k (í)+<4(í)/ k (í) = (22) 



where 



wl(t) = c\k 2 y + cj[k x - É(t)} 2 + M 2 (t) , (23) 

and E = (d/Ae)E' . We consider the case of effectively 
zero temperature, when there are no vortices in the initial 
state. So, eq. (22) is solved with the vacuum initial 
conditions 



h(t -> -oo) = [^fvr^exphiwn* - *<)] - ( 24 ) 

where V is the two dimensional volumc of the film, U is 
some initial moment of time, lü^ = [c\k 2 + M 2 ] 1 / 2 , and 
Aí = M(t-> -oo). 

Once a solution to the initial problem (22)-(24) is 
available (e.g., from a numerical integration) , one can 
obtain various quantities of interest as averages over the 
vacuum of the operators a and (3. In what follows, we 
consider three such quantities: the average vortex cur- 
rent, the energy, and the vortex occupation numbers, all 
as functions of time. 

Only the x component of the average vortex current 
is nontrivial. It can be found by averaging the operator 
expression 

A J x (t) = -i X ^d xX + i{d xX ] )x ~ 2S x t x (25) 
over the vacuum of a and /?, to obtain 



(0), 



where 



(Jx(í)> =E«( í -y . (26) 



q{t,k y ) = 2c\Y J [k x -É{t)]\h\ 2 . (27) 



If we integrate (26) over time, we will obtain the average 
vortex number transported in the x direction per unit 
length in the y direction during the entire pulse. This is 
a convenient measure of vortex transport. 

At zero temperature, and for a sufnciently adiabatic 
pulse, the only source of vortex transport is vortex tun- 
neling. In the real-time formalism, the corresponding 
contribution to the current (26) appears as a result of 
the discreteness of modes. It is exponentially suppressed 
with L x and should reproduce the result (2) of the in- 
stanton approach. 

We should be careful, however, about the ultraviolet 
regularization of eq. (27). A sharp momentum cutoff is 
not adequatc for our purposes, especially since we nced 
a regulator that would preserve the symmetry under the 
transformation (17)-(19). A suitable choice is a Pauli- 
Villars regulator — an additional field with a very large 
"mass" M', whose contribution is added to (26) with an 
opposite sign relative to that of \- If the maximal k x 
is some A 3> M', the regulator contribution to (27) can 
be computed analytically. Both the pair production and 
tunneling are negligible for large M', so we replace the 
regulator mode functions with their WKB expressions, 
and the sum over k x with an integral, to obtain 



q'(t,k y ) = ^-É. 



ttL v 



(28) 



The resulting expression for q at finite (large) value of A 
is then 



QA.(t,ky) = 2c 



A 

E 

k x = -A 



[k x - É(t)]\M 2 + q' . (29) 



and the full regularized expression is 

q(t,k y )= lim q A (t,k y ) . (30) 



IV. NUMERICAL RESULTS 

As discussed in Sect. II, it is possible to consider ge- 
ometries in which the suppression of the order parameter 
and the driving force on the vortex are entirely inde- 
pendent functions of time. For example, in the doublc- 
arm geometry of Fig. 2, the suppression is determined 
by the externally controlled current /, while the driving 
force, by the biasing flux. Accordingly, we consider here 
a situation when the order parameter is suppressed for a 
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relatively long time down to some value that allows tun- 
ncling, while the driving force exists only for a shorter 
time: its biases tunneling and leads to a nonzero value of 
the average (26). 

A convenient parametrization of the force F is ob- 
taincd by defining an effective current f related to E' 
by an effective Maxwell equation 



1 dE' 



E eF 
47r nhd 



(31) 



For a sufficiently adiabatic pulse, the first term here is 
much smaller than the second, and we neglect it in what 
follows. According to cq. (7), it is j/n s , i.e., a quantity 
akin to the vector potential, that determines how close 
the film is to criticality. So, we define an effective vector 
potential A' via the London formula 



c_A[_ 
4tt ' 



Combining (31) (with ÉwO) and (32), we obtain 

hed A' 



8« E mA 2 £ A c 



(32) 



(33) 



where «em = 1/137 is the fine-structure constant, and 
A c = $o/27rÇ is the critical vector potential. 

We consider only biasing pulses that have very small 
A' /A c ratios. Such pulses do not significantly modify A. 
So, the only difference between A and the unperturbed 
value of the penetration depth A is due to the broader 
pulse of the current I. Similarly, the vortex "mass" M(t) 
during the biasing pulse may be assumed constant and 
equal to some Mq. 

For the parameters of the film, we take d = 4 nm, 
£ = 20 nm, and A = 150 nm. The latter two vàlues take 
into account the suppression of £ and the increase in A 
due to the small value of the thickness d [14]. The sizes 
of the film are L x — 1/zm and L y = 10/xm. 

We assume that the suppression of superconductivity 
by a pulse of / has reduced n s by a factor of 25. Then, 
A 2 = 25A 2 . In our numcrical integrations, we keep A and 
the form of the pulse fixed and scan over different vàlues 
of M . 

The limiting vortex speed can be obtained by estimat- 
ing the inertial mass of the vortex m v and taking the 
ratio c 2 = hM /m v . In many cascs, the main contri- 
bution to m v comes from the small variation of electron 
density at the vortex core. In this case, m v ~ mkpd, and 
ci ~ vp [15-17]. However, when we search for an opti- 
mal tunneling path in the Euclidean time, the density 
at the core becomes a variational parameter, and it is 
advantageous for it to diffcr from the average density as 
little as possible. In this case, the main contribution to 
m v comes from the elèctric field produced by the moving 
vortex, resulting in a much larger c\ <~ (Ç/A)c [15]. For 
the above vàlues of the parameters, we use c/ci = 7.5. 



In what follows, we adopt the system of units in which 
ci = 1 and all lengths are measured in microns. Thus, 
the unit of time is 1/iím/ci = 0.025 ps. 

The easiest way to implcmcnt the quantization con- 
dition (20) is to consider A' whose integral over time is 
zero. We set 



A c t 



(34) 



and take C = 0.005 and t = 80. The latter corresponds 
to 2 ps. 

Equation (22), with 



É(t) 



cd 



8aEMA 2 £ 



L 



(35) 



and different vàlues of M(t) = Mq was integrated numeri- 
cally using a Rungc-Kutta sixth-ordcr integrator. We use 
N x -1 vàlues of k x : k x L x /2n = 0,±1, . . . , ±{N x /2 - 1), 
with N x = 32. 

Taking the limit in eq. (30) requires correcting the nu- 
merical data at least by terms of order Mq / A 2 . In our 
case, A = 307r. The correction was computed by assum- 
ing that it dominates the transport already for Mo = 10, 
which is a large enough value to significantly supprcss 
tunneling. 

In Fig. 3 we plot the total vortex transport 



Q(ky) 



1h{ky, t) + 0.0215 Mq 



(36) 



for k y = and several vàlues of Mq; tf^ — ±400. The 
data are well fit by a curve proportional to exp(— MqL x ), 
which is the instanton exponential. Note that in this 
calculation the maximal value of E is E max w 9, which 
exceeds the mode spacing 2ir/L x — 2tt. For E <§C 2n/L x , 
agreement between the instanton and real-time calcula- 
tions was confirmed analytically in ref. [10]. We now 
confirm the agreement beyond the small E limit. 
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FIG. 3. Points: total vortex transport (36) for k y = and 
different vàlues of Mq. Line: a const.x exp( — Mo) fit. 
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Because ío 3> 2it/Mq for all vàlues of Mq in Fig. 3, the 
transport is to a good accuracy adiabatic. The measure 
of adiabaticity is the vortex occupation numbers 



Tlk(í) 



V 



Wk(í) 



E k (t) - 1 



(37) 



where 



that a sufhciently slow, adiabatic variation of the biasing 
field can lead to a sizeable vortex current without any 
real vortex-antivortex pairs remaining in the final state. 
This means that induced vortex transport may be a suit- 
able technique for applications requiring a high-degree of 
quantum coherence, such as quantum memory. 



£ k (í) = |/k| 2 + ^(í)l/kï 



(38) 



are the energies (divided by K) of the individual modes. 
In Fig. 4 we plot the total occupation number 



N{k y ,t)= J2 n kW> 



(39) 



for k y = and Mo = 5. We see that there is practically 
no residual excitation (real vortex-antivortex pairs) in the 
final state: at t = 400, we obtain N(k v = 0) ~ 10" 12 . 
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FIG. 4. Total occupation number for Mo = 5 and k y = 0, 
as a function of time. 



V. CONCLUSION 
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The results of this work are two-fold. First, we have 
shown that in certain geometries (such as the doublc- 
arm geometry of Sect. II), it is possible to independently 
control the driving force acting on a vortex as it tunncls 
across a superconducting film and the suppression of su- 
perconductivity in the film (i.e., the mass of the vortex). 
Motivated by this observation, we have then considered 
a model problem of vortex tunncling in a toroidal film. 
The driving force is modeled by an effective elèctric field, 
which biases tunneling and results in a nonzero average 
vortex current. 

Second, we have used this model setup to study the 
dcpendcnce of the tunncling rate on the vortex "mass" 
(more precisely, the pair-production frequency) in the 
real-time approach, in which vortices are represented by 
a fundamental quantum field. We have confirmed the 
cxponential dcpendcnce on the "mass" found in the in- 
stanton (Euclidean) approach. We have also confirmed 
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